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ABSTRACT 


This thesis presents an analysis of a direct current 
conduction magnetohydrodynamic generator. Laminar, in- 
compressible viscous flow is considered. Experimental 
evidence is offered to show that the laminar assumption 
is valid in the presence of strong magnetic fields. 


Efficiencies are calculated for generators of various 
ieme=oe and aspect ratio, efficiency curves are presented 
and 2 method of driving the flow is considered. 


Thesis Supervisor: Herbert H. Woodson 
ioe: Prrobessor of Electrical Engineering 
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Chapter I 


INTRODUCTION 
1.1 Background 

The minimization of sound transmitted from a submarine 
through the sea to the listening devices of a potential 
enemy 1s a major problem in submarine design. A large 
fraction of the radiated noise is generated by the submar- 
ine's propulsion plant: the main engines and the screws 
and shafting 

Consequently it would by desirable to devise a propul- 
sion system having neither shaft nor screw and a minimum 
quantity of rotating machinery. Doragh* (1) has shown the 
feasibility of propelling a submarine by means of a sea 
water jet accelerated in a magnetohydrodynamic channel. 
This thesis will be a study of the closely related problem 
of providing the electrical energy necessary to drive the 
magnetohydrodynamic flow. 

The pioneer figure in the study of magnetohydrodynamic 
channel flow was Hartmann (2). In 1937 he made the first 
detailed theoretical analysis of the flow of a conducting 
fluid in a high aspect ratio channel under the influence 
of a uniform transverse magnetic field. Hartmann and 
Lazarus (3), in that same year, reported on observations 
made on mercury flow in a transverse magnetic field. To 
date the principal use of magnetohydrodynamic channel flow 
has been as electromagnetic pumps to circulate liquid 
metal coolants in atomic power installations. 

At the present time there is considerable interest in 
the development of magnetohydrodynamic power generators. 





* The bracketed numbers refer to the Bibliography which 
begins on Page 28 





Most of these proposed generators make use of a conducting 
eas or plasma as the working fluid. This type of magneto- 
hydrodynamic generating system shows promise of supplying 
iarse quantities of electric power at thermal efficiencies 
unattainable with present equipment (4,5,6,7). 


1.2 The Object and Scope of This Study 
The object of this study is the design of a liquid 


metal magnetohydrodynamic direct current conduction gener- 
ator compatible with and utilizing the same magnetic field 
as the magnetohydrodynamic pump designed by Doragh. As 

will be shown, the use of a liquid metal, instead of a con- 
ducting gas or plasma, as the working fluid permits the 
design of a complete energy conversion cycle devoid of mech- 
anical pumps and compressors. Since the basic goal of this 
design is to achieve a system free of these mechanical noise 
sources, gaseous working fluids will not be considered. 

mae channe! under consideration is rectangular in cross 
section and of high aspect ratio. The applied transverse 
magnetic field is uniform and parallel to the shorter sides 
@f the channel cross section. Furthermore, the liquid 
metal is assumed to be isotropic, homogeneous and incom- 
pressible. 

The study will be restricted to laminar viscous flow, 
since there is little or no experimental data available 
concerning turbulent flow of an electrically conducting 
mod Ima stromg transverse magnetic fleid. In fact, the 
assumption of laminar flow may be quite valid since there 
is some experimental evidence that indicates a transverse 
magnetic field delays the transition from laminar to turb- 
ulent flow. Murgatroyd (8) in his experiments with mercury 
found that laminar flow was stable when the hydraulic 
Reynolds number was less than 225 times the Hartmann number. 
Lock (9) hypothesized that laminar flow would be stable for 
a Reynolds number less than 12,000 times the Hartmann number 





of 
for larger Hartmann numbers. Since this study will be deal- 
ing only with large Hartmann numbers, both theory and exper- 
iment predict laminar flow. 





Chapter II 
PROCEDURE 


2.1 Basic Equations 

The fluid under consideration is homogeneous, isotropic, 
non-magnetizable, non-polarizable and incompressible. It is 
memeeeterized by a permeability p., an electrical conductivity 
amass density p and a dynamic viscosity n. The equations 
that govern the motion of such a fluid are the magnetohydro- 
dynamic Navier-Stokes equation. 


ppy = -srad p + V° ¥V4+IXB, (aa) 
the incompressibility condition 


diy V =O, alee) 


Maxwell's equations 


emt H= J, _ (25a 3)) 

curl & = us, (27s) 
and 

ay a = 0; (2.1.5) 
and the constitutive relation 

J =O(E + Vv X ph). (26 ) 


Where H is the magnetic field intensity, J is the conduction 
current density, E is the electric field intensity, v is the 
fluid velocity and P is the pressure. Rationalized MKS units 
are used in all equations and all equations are written with 
respect to a fixed coordinate system. 





e.2 Pressure and Velocity Relations 

The geometry of the generator to be considered is il- 
lustrated in Figure I. Fluid flows in the z-direction 
MWitn a velocity that is a function of both x and y. A 
uniform magnetic field H is applied in the x direction. 
The channel walls at y = Z w/e are highly conducting 
buses that are maintained at a potential difference by 
the external constraint imposed by the load resistance. 
The walls at x = = d/fe are non-conducting. The electric 
field E resulting from the transverse potential gradient 
is assumed to be in the negative y direction, and the 
current density J is assumed to be in the y direction. 

The channel cross section is shown in Figure II. 
Load current I is assumed to flow symmetrically in the 
external circuit. This insures that the magnetic field 
resulting from the flow of load current will be in the 
zZ direction and so will not affect the generator oper- 
ation. End effects will not be considered at this time. 

Under these assumptions, the z component of Equation 
(2.1.1) becomes 


a= 
Be, ) V. ‘ 7 
ee - ~ +,9 "2 tne V5 JB, (27 aes) 
yx sy* 


The elimination of Jy by use of Equation (2.1.6) yields 


Bay 4 yey _ Spey + © (EB Oe 0 = eye) 


a 
>x* > y< ap 1 rey 4 








This equation can be written as 





* From this point on indices indicating vector directions will 
be omaeeead. All directions are as defined in Figure f. 








Figure I 


Channel Geometry 
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Oa s <y 2 
o - kv + P=0O Cane. 
a dy ie 
where 
e C 
2 ae (d/2) _M (2.2.h) 
-(a/2)* (a/2)* 
M = B(a/2)-)= (Hartmann number) (25en5)) 
and 
p= (zp-idP ee 
m | Sos) (25270) 


ee oe , (x) cos (2n : 1)my gy (2.27) 


is desired, where 
wh 
Vv ie ‘ v cos (2n-1)ny dy (2.2. 0)) 
en-1l = ——— 
—wW L W 
If Equation (2.2.3) is multiplied by the function 


= os (2n-1)xy dy and the resulting equation integrated 


W 
W 


over the interval y = -w/2 to y = w/2, and the boundary 
condition, v(x, =w/2) = 0 is applied, the result is 








: ere oa Mo [ (enu)n ]° Von 1 + 
ot e 
Yue | (a2) , 
or 
ili 
2 2 ee) 
> Von-] - M Nee + 6) 
5 xe (2n-1)t 
where 
m2 = Me (2a hr 


Define the function 


6 = M- 


n+ii 
V ) 
cn-1l n en-l 


ee eal 
(2n-1)nx 


then Equation (2.2.10) can be expressed as 


Z 
ree of eee 2 
eal = Win Syl = 0 
yx* 
This has the solution 


Son] = Acosh M* + B sinh MX 


ee 


ek <a _ 


(2n-1)n 
(252.9) 


(252 e)) 


(2. 24258) 


(Les I 2) 


(2.2.13) 


(2ei2 in) 


where A and B are arbitrary constants to be evaluated from 


the boundary conditions. 


Because of the symmetry of the problem, v, and hence 


6 must be an even function of x. Therefore B= 0O. 


cn-l 


A may be evaluated from the requirement that v(=d/2,y) = 0. 








ime selutéion for Gon] is then 
iol oe al 
_ ~4P(-1 
oo ah WP (-1) Were. 15) 
(2n-1)n cosh M2 
and 
ier 
ae jut COREP Ee (2.2.16) 
2 
(2n-1)2M, cosh Md 
2 
whence 
vas 
a al! 
oe hp (-1) 1- cosh Mx sos (2h-1)1 ny 
Tl en a 
- (2n-1)M_© cosh M4 4 
m=] ae, 
(2ree) 


The expression for the mean velocity v can be found by 
integrating Equation (2.2.17) across the channel with the 
result 


~ 
72 ci | 1- * tanh (M |. 
— => re ee Oe =. ite 
t,o (2n-1) a Md 
(2.2.18) 


Equations (2.2.6) and (2.2.18) can now be combined to 
give the pressure gradient down the channel. 





1Se 





2 
Pe = 9° =soEp - ™ i“ 
dz = il 
[B= en 
M 
(alee pe n ‘as 
a (2n-1) M_ 


(2.2.19) 


2.3 Terminal Impedance 

For large Hartmann numbers Equation w2.2.17) predicts 
that the normal parabolic velicoty profile associated with 
laminar flow is distorted by the presence of the transverse 
magnetic field. The flow in the center of the channel is 
retarded and that near the walls acceleréted so that the 
result is a core of uniform flow and a boundary layer 
whose thickness is of the ordere of d/M. ** 

In determining the terminal characteristics of the 
generator the velocity will be assumed to be uniform and 
equal to v The Hartmann number will be large enough to 
make this assumption valid. 

With the assumption of a uniform velocity profile, 
the equivalent electrical circuit for the generator can 
be determined quite easily from Equation (2.1.6). Arter 
transforming this equation from field variables to cir- 
cuit variables the result is 


Tld . (vB + “), Ble al) 
W 


I and V are positive in the directions defined in Figure II. 





* For a complete discussion of these results see Ryabinin 
and Khozhainov (10). 
** See Shercliff (11) page 137. 





ee 


This current-voltage relationship determines the equivalent 
electrical circuit shown schematically in Figure III. 
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Figure III 


A characteristic of this generator is its low terminal 
yoltage. This is readily seen from the expression for open 
circuit voltage, vBw. For a liquid metal the maximum mean 
velocity practiable is of the order of 10 meters per second. 
The channel width is of the order of one meter. The B field 
Weed by Doragh, and hence the field acting on this generator 
is ten Webers per square meter. The open circuit voltage then 
is limited to a maximum of about 100 volts. 

2.4 Power Relations 


Define R; s the internal resistance of the generator, by 
R, = ,.. | (2 
Old 


The terminal voltage, V, can then be expressed as 


a (2.4.2) 


R, +R, 


yy . 


ame the load current, I, as 


i VvBw : (2glhes\) 


Ry ct Ry 





dye 
The electric power delivered to the load, Pi» is given by 


2 
B= Toy = (EH) FR, 


: (2.4.4) 


(R, + RL)° 


And the generator efficiency, e, is given by 





le 
e= - 
Rate of doing mechanical work on the fluid 
E C 
= 1 vBw) R 
e = - (vBw) L (25S) 
a 
Pv 1wd (R; + R,)°Pov lwd 


2.5 The Energy Conversion Cycle 


The source of energy for operating the generator is 
assumed to be a nuclear reactor. The problem of con- 
verting the thermal energy of thw reactor into mechani- 
cal energy necessary to drive the generator without the 
intervention of mechanical pumps, compressors and other 
mechanical noise sources is a difficult one. The cycle 
considered here was devised by Elliot (12). It is shown 
schematically in Figure IV. 

Elliot carries out the thermal to mechanical energy 
conversion by mixing a hot liquid metal of low vapor 
pressure with a cooler liquid metal of higher vapor 
pressure. ‘The temperature of the hot metal is above 
the boiling point of the cooler metal, and in the mixing 
process the lower temperature metal vaporizes. The mix- 
ture of vapor and liquid metal droplets is passed into a 
nozzle where the vapor expands and accelerates the liquid 
phase. The vapor and liquid phases are then separated 
from each other. The vapor is condensed and pumped back 
to the mixing chamber by means of a magnetohydrodynamic 





IES 
pump. The liquid flows through the generator, to the heat 
source and back to the mixer. 

The essential feature of this cycle is that it has no 
moving parts and should satisfy the requirement for noise- 
free operation to a high degree. Its major drawbacks are 
that it is relatively inefficient and the problem of com- 
plete separation of the two phases has not yet been 
solved. 
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Chapter IIT 


RESULTS 


3 1 Results 

Generator efficiencies were calculated for generators 
of several lengths and aspect ratios for different mean 
flow velocities. The resultant plots of efficiency versus 
velocity are shown in Figures V, VI and VII. Figures VII, 
IX and X show characteristics of a generator one meter 
wide, one meter long and with an aspect ratio of five. 

Liquidlithium was assumed to be the working fluid in 
all cases. 
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Figure VII Efficiency vs. Velocity for 
2/1 Generator 


> 





Boo 


“Joo 


50 


PmECU 


$oo 


6 xe 


Mos 


Bad |. 





19. 


\=4 
We | 
W/o<5 
Unermaun No, G.& Xo! 


Figure VIII 


Power vs. Velocity 
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Chapter IV 


DISCUSSION OF RESULTS 


fH 1 ascussion 

The results show that, with the exception of the thin 
high aspect ratio channel there was relatively little 
change in efficiency as the velocity was varied through- 
out the useful range. Also, the longer channels were 
more efficient than the short channels. This result 
was to be expected since the very low load impedence 
coupled to the generator caused large current flows. 
These currents were large enough that the viscous los- 
ses in the fluid were insignificant when compared with 
Ohmic losses. Viscous forces were important only in 
the thin channel at high speeds. 

The results also demonstrate the obvious fact that 
channel efficiency increases as channel cross section is 
increased Viscous losses are basically surface effects 
and anything that is done to decrease the ratio of sur- 
face area to volume of the channel serves to lower the 
viscous losses. 

Ohmic losses, on the other hand, are increased as 
channel width is increased and are lowered when length 
or height is increased. 

At this time it might be worthwhile to test the 
validity of the assumption of laminar flow upon which 
the whole analysis was based For liquid lithium* in 2 
magnetic fluz density of ten webers per Square meter 
the Hartmann numbers ranged from 5 x 104 to 2x 10° 
The hydraulic Reynolds number for these same flows 
ranged from 10° to 10° The ratio of Reynolds number 
to Hartmann number was at all times less than the 
Cratical valtie of 225 found by Murgatroyd. Phere is 


pps hhiee l a N ————— 


* For the properties of lithium see Appendix A 





rah 


Still seme Question as to the validity of the relationship 
for the critical Reynolds number being equal to 225 times 
the Hartmann number for the extremely high Hartmann numbers 
encountered here. In this respect it is interesting to 
mote that, in his experiments, Murgatroyd was unable to ob- 
tain turbulent flow for a Reynolds number of 10°. 

One of the major disadvantages of the constant velocity 
channel is the large pressure gradient required to produce 
any significant amount of electric energy. Since the liquid 
metal is incompressible there is no conversion of thermal 
energy to electrical energy within the channel. The kin- 
etic energy of the flow also is constant in the channel 
Thus the difference between generator inlet pressure and 
outlet pressure alone must account for both channel losses 
and power delivered to the load. 

Another disadvantage of the liquid metal channel as a 
generator is the low terminal voltage available. It is 
desirable to operate at low flow velocities to keep viscous 
losses low. Since the terminal voltage of the device is 
wucvercce by the preduct of magnetic flux density, velocity 
and channel width, a terminal voltage of no more than one 
hundred volts is expected. This low voltage requires a 
low load impedence in order to extract any significant 
amount of electrical energy from the flow. 

4.2 Conclusion 

re 1iquid metal conduction generator is extremely 
attractive because of its small size and simplicity. The 
major disadvantages listed above are serious but not in 
themselves sufficient grounds for discarding the device. 

It is only when one looks further than the generator 
itself that insurmountable difficulties appear. The 
absence of an efficient method for converting thermal 
energy to flow energy without using mechanical equip- 
ment is a severe handicap. Elliot's cycle mentioned 





eh. 
earlier has many drawbacks. It is inherently inefficient 
and its temperature requirements exceed the capabilities 
of presently available materials. The cycle requires the 
use of a liquid metal with a high boiling point in the 
liquid loop. It follows that any metal meeting this re- 
Quirement would be solid at room temperature. As a result 
a means must be provided for keeping the metal from freez- 
ing when the system is not in use. If lithium is used 
this means maintaining a temperature of above 35° EF. 

The flow could be driven by a mechanical pump opera- 
ting in a conventional steam cycle. However, this would 
defeat the major reason for having an MHD generator: 
noise elimination. Besides, if a conventional steam 
cycle must be built to operate a pump, it may just as well 
wrave a rotating eénerator directly. 





Chapter V 


CONCLUSIONS AND RECOMMENDATIONS 


5.1 Advantages 

The prime advantage of the generator studied is its small 
size and simplicity. In addition, it meets the basic require- 
ment of relatively noise-free operation. 
5.c Disadvantages 

The disadvantages of this type of generator system are 
manifold. Terminal voltage is low, losses in the system 
excluding the generator are relatively high and there is a 
serious materials problem. 
5.3.Recommendations 

Since the difficulties of producing electrical energy 
by means of a liquid metal cycle far outweigh the advantages 
of such a system, as study of an MHD power generator system 
using an ionized gas as a working fluid would certainly be in 
order. If such a system could be made small enough in phy- 
sical size to be placed in a submarine it might very well 
provide the quiet operation we are seeking. 
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APPENDIX A 
THE PROPERTIES OF LIQUID LITHIUM (13) 


Melting Point 55 ae 
Boiling Point 208° FP. 
Density .507 gm/ee 200° C. 
490 hoo 
Pera 600 
57 809 
ya 1000 
Viscosity .5918 Centerpoises 183.9° C. 
Dug 1a 2 
.55)1 208.1 
4.917 250.8 
4548 205.5 
Resistivity W5.25 micro ohm-cm 230° C. 


Hartmann Number = B(d).[o = 6.82 x 10° (a) 
a Vm 3 





APPENDIX B 
METHODS OF COMPUTATION 


ie cateutations were made using the properties of 
lithium listed in Appendix A 

Equation (2.2.19) was solved for Po, Equation (2.2.)) 
for P, and Equation (2.4.5) then solved to determine 
efficiency. 

A magnetic flux density of ten webers per square meter 


was assumed in all calculations. 
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